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We examine a rapidly rotating Bose-Einstein condensate in an anharmonic confinement and find 
that many properties such as the critical rotating frequency and phase diagram are quite different 
from those in a harmonic trap. We investigate the phase transitions by means of average- vortex- 
approximation. We find that the vortex lattice consists of a vortex array with a hole in the center 
of the cloud as the rotating frequency Q increases and the vortex becomes invisible when fi reaches 
some value. 
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a. Introduction: There have been many studies on 
Bose-Einstein Condensate (BEC) rotating with high an- 
gular momentum, of which many features resemble some 
of the novel solid-state systems such as type-II supercon- 
ductor and quantum Hall liquid. In a number of ear- 
lier experiments Q, 0, S Q where the gas is confined 
in a harmonic trap, the energetically favored state takes 
a triangular lattice of singly quantized vortices In 
the case of harmonic trapping potential, the centrifu- 
gal force would have exactly cancelled the trapping force 
and the system would have collapsed if ft reached to the 
trap frequency uj xy . Thus the centrifugal force prevents 
the rotating frequency ft from being enhanced over to xy . 
This limit makes it impossible for angular momentum 
to reach the higher magnitude as expected. To resolve 
this problem, some theoretical studies suggest that ft is 
no longer bounded by uj xy once an anharmonic term is 
introduced into the trapping potential. Vortices are suc- 
cessfully created in experiment 6| by adding a quartic 
term in the trapping potential and their phase transi- 
tions are investigate by varying the ft. The experimental 
images make the theoretical suggestion that the picture 
should become richer in an anharmonic trap more con- 
ceivable 

II SHE IH 13 Indeed, for ft > ft c (ft c 
is the critical frequency), the gas can only form an array 
of regular vortices in a harmonic trap. However, in an 
anharmonic trap, it can also be a state of multiple quan- 
tization and a mixed state except for the above state, The 
mixed state consists of a multiply quantized vortex (hole) 
at the center of the cloud and singly-quantized vortices 
around it. 

Currently, people's interest extends to the fast rotating 
BEC which exhibits richer picture. We study the fast ro- 
tating BEC in quadratic-plus-quartic trapping potential 
in this paper. First we determine the critical frequency 
ft c by means of the variational method which differs from 
that in a harmonic trap. We prove that the system in an 
anharmonic trap can form multiply-quantized state in- 
deed [l4j. Then, we extend the approach proposed by 
Ho to investigate the phase transitions between the 
mixed state and the state with an regular vortex lattices. 
Comparing the energies for different I, we obtain the ft/j 
over which a hole is generated at the center of the cloud. 
The increase of ft will make the hole to absorb the singly- 



quantized vortices around it gradually and form a giant 
vortex state finally. Because the velocity that the hole 
absorbs the singly-quantized vortices becomes larger and 
larger, the fragility of the vortex lattice increases. How- 
ever, recent theoretical study [15( found that much longer 



time were required for ft 



to reach a well-ordered 



lattice. Based on the this argument it is reasonable that 



vortex can not be detected for ft 



m experiment 



We consider a Bose-Einstein Condensate confined in 
an anharmonic trap whose potential is well approximated 
by a superposition of a quadratic and a quartic potential: 
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m 



v(f)^-Mco 2 z z 2 + -Mu; 2 r 2 + 



MuJ 2 y 4 



(A>0), (1) 



where M is the atomic mass; tu z and the frequen- 

cies of the harmonic potentials along the z-axis and in the 
xy plane r = (x,y) respectively; a±_ = [fi/(Mw I9 )]' de- 
notes the oscillator length and A a small dimensionless 
parameter characterizing the strength of the quartic po- 
tential which is A k 1CP 3 in the experiment 6J. The 
energy functional is given by 



E[V] = J d 3 rf*(h z 



hi )* 



(2) 



where \& is the condensate wave function, h z and h± refer 
to the single-particle Hamiltonian along z-axis and in the 
xy plane, namely, h z = -(h 2 /2M)V 2 z + (Mlo 2 z z 2 )/2 and 
h± = -{h 2 /2M)W\ + (Mcj 2 xy r 2 )/2 + (Mco 2 xy X^)/(2a 2 L ). 

To find the optical form of 'F for a rapidly rotating 
system, we ought to minimize the energy Eq.J2J) sub- 
ject to two constraints. Since both the total number 
of particles and the total angular momentum are con- 
stants, two Lagrange multipliers are introduced such that 
8{E — fj,N — ni z ) = 0. Here the chemical potential /i 
and the rotation frequency ft of the trap can guarantee 
the particle number as well as the angular momentum 
to be constants [l^. We therefore solve the conden- 
sate wave function by minimizing the following Gross- 
Pitaevskii functional: 



K= \ d 3 r^*[h z + h ± -riL z - 



9 - I d 3 r m 4 



(3) 



2 



in which L z is the angular momentum along the positive 
direction of z-axis ; g — (A7rh 2 a sc )/M is the strength 
of the effective two-body interaction with a sc being the 
s-wave scatting length. 

b. Critical rotational frequency- Q c : In reality, for a 
large condensate, the density profile along z satisfies the 
Thomas-Fermi approximation, which is similar to the sta- 
tionary case. Since the gas rotates around the z-axis and 
its feature in the xy plane is more interesting, we simply 
assume that the particle density along z is a constant p 
|16| so that we have a two dimensional quantum system 
(e.g., z = plane without loss of generality). 

It is instructive to consider the limit case g — and 
A = when the hamiltonian in xy plane is similar to that 
in the quantum Hall regime. In this case the eigenen- 
ergy gives rise to the Landau Levels for which the wave 
function <I> m for the lowest Landau Level (LLL) reads: 
$ m oc r m exp(imcj)) exp(r 2 /2a^_) ^?|]- in experiment, 
the g and A are nonzero but their magnitudes are rela- 
tively small, thus we can solve the problem by means of 
the variational method. Let us write the variational wave 
function as <I> = Cr m exp(im^) exp(P 2 r 2 /2a 2 L ) with vari- 
ational parameter (3. Here C = (3 m+1 /[a , l(jia 1 2 m\) 1 / 2 ] 
is the normalization constant. The energy for state 
^ = yfpQ is given by (in units of phu> xy ) 



K 



<> 2 



P 2 ^ 



)ri(m) 



n 



_A 

2/3 



4 r 2<» 



+pa sc (3 r 3 (m) + ( p), 



(4) 



where Ti(m) = f3 2 (r 2 ) ja\ = to + 1 refers to the rms 
width of the wave function in units of P 2 /a 2 Ll ^(m) = 
/3 4 (r 4 )/a 4 L — (m + l)(m + 2) the expectation value of 
r 4 , and T 3 (m) = 2^a\/ p 2 / d 2 r | $ | 4 = 2m!/(2 2m (m!) 2 ) K(z) 
the interaction contribution. Minimizing the energy with 
respect to the variational parameter (3, we obtain the 
following two equations with respect to j3 for the cases of 
to = (non-vortex configuration) and to = 1 (one vortex 
of unit strength), respectively, 



P$ - p 2 - 4A + 2pa sc p 6 
2p\ - 2p\ - 12A + pa sc {3\ 

Their corresponding energies are 



0, 
0. 



(5) 



K o = + i) - £-) + A x 2 + pa scf 3 2 

+( m), 

K i = (^(/3i 2 + i)-— )x 2 + ^4x6 + P a sc /3 2 



Pi' 



xy 



2pf 



+ ( M)- 



(0) 



The experiment data N = 3 x 10 5 , aj xy = 2tt x 65.6Hz, 
2 = 27r x 11.0Hz, and a sc = 53 A, give rise to pa sc w 20 



[T^| . Thus we have approximately po = 0.397664 and 
/3i = 0.551786. For K Q =K X , we obtain fl c /uj xy = 0.2223 
with O c the critical frequency. The state without vor- 
tex is stable when f2 < f2 c , but when £1 > f2 c the state 
with vortex is stable. For A = 0, the critical frequency 
is Vt c /u) xy — 0.136 which is smaller than that for A ^ 0. 
It has been shown 18] that the gas without interaction 
in a harmonic trap does not generate vortex. The con- 
tributions of both interaction and rotating frequency f2 
result in vortex in the condensate. The larger the value 
A is, the smaller the contribution of the interaction term 
in the total energy will be. Thus, larger rotating fre- 
quency SI is required for generating vortex. We examine 
further state configurations with larger units of vortex, 
viz., for to = 2, 3 • • • [T^ . The energy K2 is minimized 
when Po = 0.642021. K Y = K 2 gives fl 2 /uj xy = 0.7128. 
Clearly, when fl2/u>xy > 0.7128 the K\ becomes larger 
than K 2 , which implies that the state with doubly quan- 
tized vortex is more stable than the singly quantized vor- 
tex. Unlike the case in harmonic trap where the state 
with singly quantized vortex (unit quantum number) is 
merely most favorable when ft > il c , the stable state 
in a anharmonic trap favors vortex with large quantum 
numbers (angular momentum). The faster (the larger 
Q) the anharmonic trap rotates, the higher the quantum 
number of the vortex will be. 

c. Vortex lattice: We consider the state with a large 
angular momentum by looking at configurations with a 
large number of vortices of unit strength centered with 
a giant vortex of larger strength. Since the dynamics 
along z-axis is identical to the stationary case, the kinetic 
energy | V z $ | 2 in z direction is neglectable in Thomas- 
Fermi approximation (TFA). Let us write K = J dzK(z), 



d 2 r^* 



2to 
1 



vi 



^A, 4 
2a 2 ; 



Q.L, 



f d 2 r\^ 



(7) 



where p(z) = fi — Mw 2 z 2 /2. We are now in the position 
to calculate K(z). The condensate wave function is ex- 
pressed as $ = J2m=o Cm^m ^ ex p( — r 2 /2a 2 ± )F(u) with 
F(u) = Ttf =1 [u — b a ], the vortex function of u. Here u 
refers to (x + iy)/a± and b a the position of vortex 
A hole in the center of the gas was shown to be gen- 
erated with enhancing f2 in weak repulsive interaction 
regime 0, which was also notice in numerical simula- 
tion p| . Hence we assume that {b a } form a regular lattice 
in which the vortex at the center site carries I units of 
angular momentum while the vortices at the other sites 
merely carry a unit of angular momentum. Then F(u) 
can be written as 



F(u)=u 1 - 1 l[(u-b a ), 



(8) 



where 6q = is implied. The general form of ^ is likely 



* = /(*)*, 



d 2 r 



3| a =l, 



(9) 



3 



where f(z) describes the density profile of the gas along 
z-axis and <& = <t/D(z) with 



<i> = cxp(-r 2 /2ai)u l - 1 Y[(u~b a ), 



(10) 



1=1 



and D 2 



J d 2 r | $ | 2 , the normalization constant for 



Now the number constraint J d 2 r \ ^ \ 2 becomes 
/ dz | / | 2 = iV. Noting that in the complex coordinate 

L z = h(udu — u*du*), 

we can obtain the expectation value of L z , 



**L Z * = h 
Then Eq.Q becomes, 



K(z) = 



+ 



-p,(z) - 



[(r/a ± f 1] | * | 
h(uj xy - VL)(r 2 )$ 



(11) 



2a 2 



(12) 



where /i(z) = ^-ffl-(Mw^ 2 /2), (r 2 )$ = / d 2 rr 2 | 3> | 2 , 
(r 4 )<£, = )d 2 rr 4 | $ | 2 and J # = /d 2 r | $ | 4 . The 
minimization of the energy is performed by the variation 
of the parameter {b a }. 

To evaluate (?" 2 )*, (?~ 4 )* and 7$, we note that, 



$ | 2 = e 



H = H 1 +H 2 



(13) 



where = r 2 /a 2 L — 2 J^ Q ln | r — 6 Q | and if 2 = 

— 2(Z — 1) In | r | . By means of the Fourier transform, one 
obtains 0, 



Viifi = -z > cosfc-r, 



$ r = 



fc#0 

r 2(i-i) e -(r 2 /^ 2 ) TJ exp(C fe cos£r), 



7T 

■ 



4tt 



(14) 



In the average- vortex approximation the |<I>| 2 be- 
comes |$| 2 = r 2 ^ 1 ' cxp — (r 2 /cr 2 ) = /i(r) exp — (r 2 /cr 2 ); 
here /i(r) = r 2 '^ 1 ) is a slowly varying envelope func- 
tion, which gives rise to the components in higher 
Landau Levels |2l|; the $ then becomes | $ \ 2 — 
( r 2(/-i) exp (_ r 2/ CT 2jj ^/ ( 7r r(/) C r 2Z ). Using the above re- 
sult, we can easily evaluate the values. (^ 2 )$ = 



k 



1(1 + l)o- 4 , and h 



(2i-l)!/(2 2, - a [(Z-l)!] 
into Ea l|12l) . we have 

K(z) = -fl(z) + 



Substituting these results 



lh(uj X y — fl)a 2 



2a 4 



-Z(Z + l^ 4 / 2 + 



Minimizing the energy, we obtain two equations, 



dK{z) 



= 2 



-fl(z) + 



ih(uj xy — n)a 2 



Of 

dK(z) 2lh(cu xy — Q)a „ 2hu) x „X 



^a + iy 4 / = o, 



9<T 



3 j- 2 



2hw xy a 



re = o. 



(16) 



In the light that the above nonlinear equations are dif- 
ficult to solve analytically, we study those equations for 
the cases of A = and A ^ so as to capture some 
physical implications. 
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FIG. 1: (color on line) The phase diagram in the plane of 
jjb(z) versus (1 — Q./uo xy ) which is divided by the parabolic 
curve fi(z) — —9hco X yl(l — Q/u> xy ) 2 /10(l + 1)A, the oblique 
line jl(z) = huixy [2.5Z(Z + l)A + 3/(l — Q/aj xy )\ and the vertical 
line p,(z) = —5hco xy (l + l)A/3. 

For A = 0, viz., the gas is confined in a harmonic trap, 



o 3lh(uj xy — f2) 



K/ 2 ) = 



9l 2 H 2 U! X y(uJ X y - fi)£' 

Substituting Ea. i|TT|) into Eq.JI3}, we obtain, 



(17) 



[*(*)]„ = 



[a(z)] : 



/ 2_ J_\ 

27h 2 aj xy (u! X y - Q)a sc V l 2 £ P£J 



(18) 

in which a = (— p^ + jlj) for example, [cv]; = i = —1, 
[a]i=2 = -3/4, [a] l=3 = -40/81 etc. This implies 
that the energy of the system becomes larger with in- 
crease of I. Thus the hole in the center of the sys- 
tem can not be formed in a harmonic trap, which 
also be confirmed by variational wave function $ = 
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FIG. 2: (color on line) The curves of K versus Q/u) my in unit 
of (2hu> xy ) 3 ^ 2 M -1 / 2 j (2uj z a 3C ) for different i with parameters 
2tt x 65.6, uj z = 2tt x 11.0, a sc = 53A. 



Cr m exp(«TO(/)) exp(/3 2 r 2 /2a^) by regarding (3 as varia- 
tion parameter. 

For A ^ 0, viz., the gas is confined in an anharmonic 
trap, 



a i / ± 



3(1 - Q/L0 X y) 

5(Z + 1)A 



± 



^9(1 - n/to xy ) 2 + [10(1 + l)A/i(z)] /(lhu xy ) 



5(Z + 1)A 
1(1 - Q/uj xy )a 4 A/(7 + l)cr 6 



+ 



(19) 



Since a > a±, the condition for the existence of solu- 
tion of the above equations is plotted in FigQ] In the 
regime I and II, only the root of (a 2 /a 2 L ) + satisfies the 
condition a 2 > a\. In the regime III, both (a 2 /a 2 L )± sat- 
isfy that condition, which implies the system has double 



stable states and two possible configurations of vortex 
lattice can be formed. The vortices in this regime be- 
come invisible for it is easy to vary between those two 
configurations due to fluctuation. In the other regimes, 
(<r 2 /a 2 L )± do not satisfy that condition. In Thomas- 
Fermi approximation, /i = MR 2 u> 2 /2, Ri stand for the 
maximum extents of condensate and for the frequen- 
cies in three direction |18| . According to the experimen- 
tal data [101, the strait line Jx(z) = fi — Ml — Muj 2 z z 2 /2 
must cross the regimes I, II and VI but not cross the 
other regimes. Thus the condition a 2 > a\ becomes 
fi(z) > fkJ xy [2M(l + 1)A + 3/(1 - £l/w X y)]. Note that 
for the same rotating frequency f2 with the growing of I, 
the maximum extent of condensate along z-axis R z will 
decrease, 

R 2 , oc n/hoj xy -n/uj xy -2.5l(l+l)\-'3l(l-n/uj xy ). (20) 

For simplicity, we assume straight vortex lines \J\ so that 
the (2D) density of vortex remains a constant along z- 
axis, then we obtain K = 2R Z K(Q). 

We plot the values of K versus D,/uj xy for different I 
in Fig|21 We can see that with increasing Q, the state 
becomes a mixed phase with a multiply quantized vor- 
tex (hole) at the center of the cloud and singly-quantized 
vortices around it 0|- This assures that a fast rotating 
Bose-Einstein condensate confined in a quadratic-plus- 
quartic potential can generate a hole in the center of the 
gas which is in agreement with the theoretical and exper- 
imental studies ((HH- Our calculation shows that if Q, in- 
creases slowly, an array of vortices is formed at first (just 
like in harmonic trap); then a hole carrying two unit of 
angular momentums is generated when Qh/u X y ~ 0.716, 
which keeps until it reaches fl/u xy ~ 1.084. With the in- 
creasing of f2 the hole will absorb the singly-quantized 
vortices around it whose angular momentum becomes 
larger and larger. Looking at the positions of intersec- 
tions in the three graphs in FigEl it is easy to find that 
the velocity of absorbing the singly-quantized vortices 
becomes larger and larger, and finally a giant vortex car- 
rying all the angular momentums will be formed. When 
Q/u xy > 1.084, every configurations of the vortex lattice 
remains for a short time, so it is difficult to detect the 
vortex. Our result can also explain why no vortex lines 
are detectable in the regime when reaches to some large 
value @. 

d. Conclusion: Employing the approach of Ref. , 
we obtained the new form of the condensate wave func- 
tion which is essentially the Lowest Landau Level (LLL) 
wave function with a regular lattice of vortices multiplied 
by a slowly varying envelope function. Using the varia- 
tional trial function, we obtained the critical value tth 
and proved that the fragility of the vortex lattice could 
be enhanced by increasing f2. This result is in agreement 
with the experiment result [2(j where the gas was trapped 
in a harmonic potential. In the experiment po|. the vor- 
tex lattice remains ordered, but its elastic shear strength 
is drastically reduced by increasing fl. This implies that 
even a minor perturbation to the cloud can cause the lat- 



5 



tice to melt. In the regime with very small A, the quartic 
term also provide a perturbation to destroy the vortex 
lattice, it is therefore invisible when the f2 reaches some 



large value @|. 
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